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0.1  Abstract  of  Research  Supported. 

For  the  most  part,  the  work  of  this  project  has  concentrated  on  non-convex  minimiza¬ 
tion  problems  in  nonlinear  elasticity  theory,  and  the  existence  of  deformation-induced 
microstructure  under  severe  loading  conditions.  The  following  papers  have  been  com¬ 
pleted  with  the  acknowledgement  of  the  grant  DAAL03-91-G-0184: 

[1]  Fosdick,  R.L.  t  Y.  Zhang.  The  torsion  problem  for  a  nonconvex  stored  energy 

function.  Accepted  for  publication.  Anh.  for  Rational  Mech.  Anal.(l992). 

[2]  Dunn,  J.E.  Ic  Roger  Fosdick.  The  Weierstrass  condition  for  a  special  cltiss  of 

elastic  materials.  Accepted  for  publication.  J.  Elasticity{\992). 

[3]  Fosdick,  R.L.  k  Y.  Zhang.  Coexistent  phase  mixtures  in  the  anti-plane  shear  of 

an  elastic  tube.  Submitted  for  publication.  ZAAfP{\99Z). 

In  addition,  some  resources  of  this  grant  were  used  to  help  support  the  Ph.D.  can¬ 
didate  J.H.  Yu,  whom  I  have  advised,  during  the  course  of  his  thesis  research.  This 
work  is  not  yet  completed  but  it  will  be  finished  by  June  1993,  and  Mr.  Yu  is  expected 
to  be  awarded  a  Ph.D.  degree.  The  thesis  title  is: 

[4]  Yu,  Jang-Hurng.  The  Nonlinear  Oscillation  of  Viscoelastic  Cylindrical  and  Sphe¬ 

rical  Shells,  1993. 
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0.2  Report  on  [1]. 

Motivation  for  this  work  is  contained  in  the  interesting  fundamental  experimental 
observations  of  Bader  and  Nadai(1927),  Takenaka(1928),  and  Nakanishi(1928)  con¬ 
cerning  the  torsional  deformation  of  steel  bars.  Some  pictures  of  these  observations, 
taken  from  a  classic  book  of  Nadai(1950)  on  the  flow  and  fracture  of  solids  are  con¬ 
tained  in  the  figures  at  the  end  of  this  report.  These  figures  show  that  as  the  angle 
of  twist  is  increased,  the  internal  deformation  pattern  in  a  steel  bar  eventually  un¬ 
dergoes  a  transition  in  which  many  small  sliver-like  zones  of  large  shear  occur.  These 
zones  are  generally  sector-shaped  with  a  periodic  circumferential  distribution  when 
viewed  at  in  a  right  cross  section.  Sometimes  the  zones  are  seen  to  be  wedge-shaped 
and  distributed  in  the  axial  direction,  in  a  longitudinal  cross  section  of  the  bar.  The 
experimental  observations  show  that  with  increasing  angle  of  twist  the  distribution 
of  zones  of  large  shear  becomes  more  dense  and  that  eventually  the  cross  section  of 
the  bar  contains  essentially  three  regions  of  shear  behavior:  an  inner  core  of  classical 
torsion  with  small  shear,  an  outer  annular  region  of  relatively  large  shear,  and  a  mix¬ 
ture  region  in  between,  in  which  thin  sector-shaped,  or  wedge-shaped  zones  of  large 
shear  are  distributed  periodically  among  zones  of  small  shear. 

In  [1],  we  studied  this  phenomenon  within  the  framework  of  minimization  theory 
in  nonlinear  elasticity  theory  for  a  non-convex  stored  energy  function,  and  in  so 
doing  we  developed  a  new  perspective  on  its  origin.  From  an  energetic  point  of 
view,  the  optimal  deformed  configuration  of  an  elastic  body  is  one  which  associates 
with  the  lowest  total  potential  energy  among  all  admissible  states.  Thus,  if  one 
considers  the  experimentally  observed  sector-shaped  or  the  wedge-shaped  zones  of 
large  shear  as  being  of  fundamental  importance  in  establishing  the  optimal  deformed 
configuration  of  highly  twisted  cylinders,  then  our  goal  was  to  provide  a  theoretical 
basis  for  why  there  exists  a  minimizing  configuration  for  the  total  potential  energy 
with  such  complicated  microstructures. 

Because  the  existence  of  zones  of  large  shear  distributed  among  the  zones  of  small 
shear  is  evidence  that  there  are  two  phases  -  large  shear  and  small  shear  phases  - 
in  equilibrium  coincidcntly,  we  adopted  the  idea  that  the  stored  energy  function  has 
a  loss  of  local  convexity  in  parts  of  its  domain  of  definition.  Then,  in  particular,  we 
studied  the  non-convex  minimization  problem  associated  with  the  torsion  of  a  solid 
incompressible,  isotropic,  and  homogeneous  elastic  bar  of  fixed  height  and  right  circu¬ 
lar  cross  section.  We  found  that  for  an  appropriate  class  of  admissible  deformations  in 
the  Sobolev  space  where  p  >  2,  the  classical  equilibrium  deformation  field  is  a 
minimizer  of  the  potential  energy  when  the  applied  angle  of  twist  is  either  sufficiently 
small  or  sufficiently  large.  However,  for  a  broad  range  of  moderate  twists  we  found 
that  there  is  no  minimizer,  and  for  these  cases  we  construct  a  minimizing  sequence 
which  converges  in  a  weak  measure-theoretic  sense  to  the  infimum  of  the  potential 
energy.  This  limiting  sequence  exhibits  an  infinitely  fine  microstructure  analogous 
to  a  two-phase  mixture  of  large  and  small  states  of  shear.  The  volume  fraction  of 
large  to  small  shear  is  found  to  depend  on  the  radial  position  in  the  cross-section. 
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with  the  greater  fraction  of  large  shear  being  at  the  outer  radius  and  decreasing  with 
the  radial  position  to  zero  at  a  particular  radius  determined  by  the  material  and  the 
overall  angle  of  twist.  The  general  picture  of  our  construction  is  much  like  the  above 
experimental  observations;  an  annular  ring  caries  a  dense  radial  dependent  mixture  of 
two-phase  shears  which  changes  continuously,  at  a  determined  radius,  into  a  classical 
small  shear,  single-phase  inner  core.  The  radius  of  the  inner  core  decreases  as  the 
angle  of  twist  increases.  We  found  that  as  the  microstructure  develops  and  becomes 
more  dense  with  increasing  angle  of  twist,  the  material  is  weakened  and  the  torsional 
moment  response  becomes  softer. 


0.3  Report  on  [2]. 

A  nonlinear  hyperelastic  material  is  characterized  by  a  stored  energy  function  W(*)  : 
P  —»  iR,  where  P  is  an  open  subset  of  the  set  of  second  order  tensors  on  with 
positive  determinant.  It  is  required  that  H'(-)  be  framt  indifferent  or  objective,  i.e., 

W(QF)  =  W(F)  (1) 

for  all  F  G  P  and  for  all  proper,  orthogonal  tensors  Q.  Since  (1)  requires  that 
QP  C  P,  objectivity  is  a  restriction  on  both  W(')  and  its  domain  P.  We  assume 
that  VV(*)  is  twice  continuously  differentiable  on  P. 

Let  n  be  an  open,  bounded,  and  connected  region  of  JR®  with  boundary  dfl. 
Let  y’(-)  :  d(l  -*  be  pre-assigned  and  consider  the  set  9  of  one-to-one  functions 
y(  ) :  n  — »  R®  which  are  continuous,  piecewise  smooth,  meet  Vy(f^)  Q  ^^ve  a 

continuous  extension  which  agrees  with  the  boundary  data  y*(-)  on  dfl.  In  nonlinear 
elasticity  theory  a  fundamental  role  is  played  by  functions  y(-)  €  9  which 

minimize  j  W(vy)(fu  on  O'.  (2) 

n 

As  is  well-known,  wherever  such  minimizersy(-)  arc  sufficiently  smooth,  it  is  necessary 
that  the  WeicrsJross  condition 

W(vy(x))  +  {F-  vy(x))  •  IVj.(vy(x))  <  IV(F),  (3) 

be  satisfied  for  all  F  €  P  such  that  F  —  Vy(x)  is  rank  one. 

Let  us  say  that  the  function  W{-)\'D  -*  ft  is  rank  I  convex  atF  if 

W{F)  -I-  (F  -  F) .  Wj.(F)  <  W(F)  (4) 

for  all  F  such  that  F  —  F  is  rank  t.  Also,  let  Zy*  C  P  denote  the  set  of  tensors 
at  which  W(-)  is  rank  1  convex.  From  (3),  we  see  that  the  gradient  of  a  minimizing 
field  y(-)  must  meet  Vy(^^)  ^  because  of  this,  it  is  common  to  expect 

non-smooth  minimizing  fields  whenever  1^*  is  not  comn  cted. 
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With  few  exceptions,  the  implications  of  rank  1  convexity  and  the  structure  of 
remain  important  open  problems  that  require  working  through.  In  [2],  we  present 
necessary  and  sufficient  conditions  for  W(*)  to  be  rank  1  convex  at  a  point  F  €  I?  for 
the  special  class  of  isotropic  elastic  materials  for  which 

W{F)  =  0(1),  (5) 

where  I  =  |Fp.  We  assume  that  0(*)  :I  R\s  twice  continuously  differentiable  on 
J={I|I  =  |F|^F€P}. 

As  a  rndn  result  of  this  work,  we  prove  the  following  theorem  and  study  its 
consequences: 

Theorem.  Suppose  (5)  holds.  For  F  6  P,  sef  B  =  FF^,  I  =  tr(B)  =  |Fp,  and 
let  $maz  denote  the  maximum  eigenvalue  of  B.  Then,  the  stored  energy  function 
W{-) :  V  IR  is  rank  1  convex  atF  £T>  if  and  only  if 

0'(i)  >  0, 

and  _ 

0(i)  +  K/i  +  (i-i)/^„ 

for  a// 1  €  J  such  that  I  >  i  — 

Although  the  materials  that  we  analyze  here  are  too  special  to  be  of  general  applica¬ 
bility,  they  reveal  clearly  that  subtle  but  important  consequences  flow  from  both  the 
structure  of  V  and  the  structure  of  P’’*.  In  this  work,  we  are  motivated  by  elasticity 
theory,  stability,  and  coexistent  phase  structures,  and  because  of  that  we  admit  only 
deformations  and  their  gradients  that  preserve  orientation.  When  this  essential  phy¬ 
sical  requirement  is  not  enforced  the  results  can  be  fundamentally  different  from  those 
of  (6),  as  is  evidenced  by  the  conclusions  of  Theorem  1.10  in  Dacorogna’s  impressive 
1989  book  on  the  calculus  of  variations. 


(6)i 


-i)0'(i)< 
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0.4  Report  on  [3]. 

In  (3],  we  have  studied  the  regularity  properties  of  a  solution  to  an  anti-plane  shear 
non-convex  minimization  problem  and  determined  the  existence  and  general  nature 
of  those  sub-regions  of  the  body  (bordered  by  free  boundaries  which  must  be  found 
as  part  of  the  solution)  in  which  the  strain  field  is  expected  to  be  highly  irregular. 
In  general,  the  problem  has  no  minimizer  but  we  construct  a  weakly  convergent 
minimizing  sequence  of  deformations  in  a  certain  .Sobolev  space  for  which  the  total 
potential  energy  converges  to  its  infimum.  We  show  that  the  minimizing  sequence 
has  a  weak  convergence  property  in  a  measure-theoretic  sense,  and  we  characterize 
the  regularity  properties  of  this  sequence. 
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More  specifically,  we  consider  an  isotropic  elastic  tube  whose  cross  section  is  a 
convex  ring  whose  outer  lateral  boundary  is  fixed  and  whose  inner  lateral  surface 
is  given  a  uniform  displacement  along  its  axial  direction.  Plasticity  theory  is  com¬ 
monly  employed  in  the  study  of  problems  that  involve  the  damage  of  materials  and 
the  localization  of  deformation  (i.e.,  shear  bands).  Here,  however,  we  provide  a  new 
perspective  on  the  understanding  of  the  mechanisms  of  material  damage  and  the  lo¬ 
calization  of  severe  deformations  from  an  energetic  (minimization)  point  of  view.  The 
non-existence  of  a  minimizer  to  our  characterizing  non-convex  minimization  problem 
for  a  certain  range  of  prescribed  axial  displacement  of  the  inner  lateral  surface  implies 
that  among  all  ‘‘admissible”  deformations  there  is  none  for  which  the  values  of  the 
stored  energy  function  correspond  to  its  convex  points  almost  everywhere  in  the  body. 
Because  of  this,  we  find  that  to  reach  the  inhmum  the  tube  divides  into  three  cha¬ 
racteristic  subdomains-one  of  high  strain,  one  of  low  strain,  and  one  of  intermediate 
“mixed”  strain.  In  the  intermediate  “mixed”  strain  subdomain,  the  stored  energy 
function  looses  its  strict  convexity.  The  main  variational  problem  then  gives  rise  to 
a  free  boundary  value  problem  in  which  the  subdomain  where  the  strict  convexity  of 
the  stored  energy  function  breaks  down  must  be  determined  as  part  of  the  solution. 
The  description  of  this  intermediate  phase  mixture  region  was  one  of  the  goals  of  this 
work.  Generally,  we  find  that  it  is  a  thin  localized  region  through  which  the  shear 
strain  suffers  a  large  change.  This  suggests  the  idea  of  a  shear  hand-&n  idea  more 
commonly  encountered  in  studies  involving  classical  plasticity  theory. 

As  a  complement  to  our  theoretical  development,  we  implemented  an  FEM  code, 
developed  originally  for  the  Cray  with  a  direct  solver,  and  revised  it  so  as  to  be  com¬ 
patible  with  the  Connection  Machine  parallel  environment  and  a  GMRES  iteration 
scheme.  This  code  has  a  large  refinement  capability  which  allows  us  to  study  those 
regions  of  the  body  where  the  strain  field  is  rapidly  varying  and  becoming  highly  irre¬ 
gular,  i.e.,  the  part  of  the  body  which  supports  a  fine  phase  mixture.  Our  numerical 
results  on  the  CM-200  have  established  some  interesting  global  features  of  the  shape 
of  the  boundary  of  these  important  regions  (free  boundaries)  that  have  so  far  alluded 
mathematical  study.  They  also  have  confirmed  much  of  our  theoretical  analysis  of 
this  problem.  The  fine  phase  mixture  regions  have  the  appearance  of  shear  bands, 
and  they  are  responsible  for  the  dramatic  softening  response  of  the  body  as  the  loads 
are  increased.  Several  additional  exploratory  numerical  results  are  reported  in  [3] 
for  this  problem  concerning  bodies  of  various  geometric  cross  sections  and  materials 
with  stored  energy  functions  having  different  nonconvexities.  As  an  extension  of  both 
[1]  and  [3],  we  have  begun  to  investigate  the  consequences  of  incorporating  regulari¬ 
zing  higher  strain  gradient  terms  into  the  constitutive  equation  so  as  to  penalize  the 
emergence  of  phase  mixture  zones  and  transition  layers  of  rapid  strain  variation. 
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0.5  Report  on  [4]. 

In  this  thesis,  Mr.  Yu  has  studied  the  relationship  between  differential  type  and 
history  type  constitutive  equations  for  describing  the  nonlinear  behavior  of  solids. 
While  the  fading  memory  approximation  to  history  functionals  for  retarded  motions 
is  known  to  lead  to  differential  type  constitutive  relations,  often  this  is  assumed  to 
imply  that  the  solution  of  equivalent  problems  for  the  two  constitutive  theories  will 
be  “close”  to  one  another  in  some  reasonable  sense.  But,  this  is  not  generally  true. 
A  part  of  this  dissertation  is  concerned  with  this  issue  within  the  context  of  the  non¬ 
linear  motion  of  solid  viscoelastic  cylindrical  and  spherical  shells.  The  fundamental 
balance  laws  for  this  work  are  of  two  kinds-second  order  nonlinear  ordinary  diffe¬ 
rential  equations,  and  second  order  nonlinear  ordinary  functional  integral-differential 
equations  of  the  delay  type  (highly  non-standard).  The  elastic  contribution  to  these 
equations  is  generally  due  to  a  non-convex  strain  energy  function  which  carries  with 
it  the  possibility  of  multiple  equilibrium  states,  and  this  considerably  complicates  the 
story.  In  this  case,  chaotic-like  motions  have  been  determined. 

The  fading  memory  approximation  for  retarded  motions  has  played  an  important 
part  in  the  development  of  constitutive  theory  in  continuum  mechanics  since  the  early 
1960*8,  but  it’s  practical  consequences  are  far  from  understood.  In  fact,  since  only 
constitutive  theory  is  invoked  in  this  approximation,  it  is  not  clear  that  adopting  this 
procedure  is  at  all  reasonable  if  approximate  solutions  of  the  equations  of  motion  are 
what  one  is  after.  A  study  such  as  that  undertaken  in  [4]  has  a  significant  bearing  on 
the  development  of  approximation  theory  for  the  solution  of  problems  in  non-linear 
continuum  mechanics. 
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Figure  1:  Cross  sections  of  soft-iron  bars  of  circular  cross  section  subjected 
to  a  torsion  deformation.* 

J.O.OOT  ».  }»0.04  *,/«,  «.  *»/«. 


Figure  2:  Cross  sections  of  soft-iron  bars  with  a  square  cross  section  sub¬ 
jected  to  a  torsion  deformation. 

jJ=OJS  ^  ^=0U4t,lcm.  c.  jj=0.90  *,/*-!. 
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Fii^ure  3:  Cross  sections  of  soft-iron  bars  with  a  rectangular  cross  section 
subjected  to  a  torsion  deformation. 


Figure  4:  Cross  sections  soft-iron  bars  with  a  triangular  cross  section 
subjected  to  a  torsion  deformation. 
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